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Abstract Bubble-driven inertial pumps are a novel m- 
ethod of moving liquids through microchannels. We co¬ 
mbine high-speed imaging, computational fluid dynam¬ 
ics (CFD) simulations and an effective one-dimensional 
model to study the fundamentals of inertial pumping. 
For the first time, single-pulse transient flow through U- 
shaped microchannels is imaged over the entire pump 
cycle with 4-/is temporal resolution. Observations con¬ 
firm the fundamental N-shape flow profile predicted 
earlier by theory and simulations. Experimental flow 
rates are used to calibrate the CFD and one-dimensio¬ 
nal models to extract an effective bubble strength. Then 
the frequency dependence of inertial pumping is stud¬ 
ied both experimentally and numerically. The pump ef¬ 
ficiency is found to gradually decrease once the succes¬ 
sive pulses start to overlap in time. 

Keywords Micropumps • Inertial pumps • High-speed 
imaging • CFD • ID models 

PACS 47.60.Dx • 47.61.Jd 

1 Introduction 

Continuing development of microfluidics requires reli¬ 
able integrated components such as micropumps, mi¬ 
crovalves, micromixers and microheaters P3[2j[3llll[5] . Of 
particular interest are devices with no moving parts 
because of their relative robustness and ease of fab¬ 
rication. Such dynamic devices have to operate under 
different principles than their more conventional me¬ 
chanical counterparts. One approach has been to utilize 
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high-pressure vapor bubbles as mechanical actuators. 
The bubbles can be created by compact microheaters 
embedded inside microchannels and they have enough 
power to displace fluid in channels as narrow as sev¬ 
eral microns across. The challenge then is to convert 
displacements into useful microfluidic functions. 

Several bubble-driven micropump designs have been 
proposed. Some of them utilize asymmetric nozzle-dif¬ 
fuser geometries ISlITl while others use multiple heaters 
to create asymmetric heating or a traveling vapor 
plug [iniiiiiiniiniiTi]. a more recent development is the 
inertial pump that relies on dynamic asymmetry of the 
bubble’s expansion-collapse cycle when a microheater 
is located close to one end of a channel [TKlITHlITTinin 
mini. Such pumps have characteristic sizes from sev¬ 
eral microns to hundreds of microns. They can operate 
even in straight rigid channels of constant cross-section 
made of robust materials. These highly integrated de¬ 
vices can be made by the batch fabrication processes 
used in large-scale production of inkjet printheads [21]. 

We first reported fabrication and characterization 
of highly integrated inertial pumps in Ref. [T7|. This 
was followed by a detailed theoretical analysis of the 
pumping principle in Ref. m- The models predicted 
unusual N-shape flows within 50-/is long pump pulses. 
In the first experimental paper m this feature was 
not reported. Resolving it experimentally required high¬ 
speed cameras with frame rates of up to a million frames 
per second and micron-scale spatial resolution. This is 
the subject of the present work. 

Even without intricate details of bubble generation 
and interaction in solid-wall tubes [2011^123] . the phy¬ 
sics of inertial pumping is complex. Net flow is a result 
of a subtle imbalance of mechanical momenta of two 
fluidic columns colliding when the vapor bubble col¬ 
lapses HEIIITIIIHI- The imbalance is caused by unequal 
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Fig. 1 (a) A high speed camera image of tested devices with a vapor bubble near maximum excursion. A dark black area in 
the bottom of the micrograph is a fluidic reservoir, and a grey boundary near it is a “shelf” of the channels. Channel width 
is 22 /im, channel height (out of the image plane) is 17 /rm. (b) A corresponding plan view. The dark rectangle is a pump 
microheater. The scale is she same as in panel (a). 


fluid inertia during bubble expansion and asymmetry of 
vortex formation near channel-reservoir interfaces m- 
Fluid in the shorter channel arm has lower inertia and 
reverses its direction earlier. It has more time to accel¬ 
erate under the external pressure and acquires a larger 
mechanical momentum by the time of collapse. At col¬ 
lapse, the momentum excess is converted into uniform 
flow of the entire channel. The latter proceeds in the 
direction from the short channel arm toward the long 
one until dissipated by viscous forces. In well-designed 
pumps, displaced volume per pulse can be as high as 
40% of the maximal bubble volume. At ^ 10^ Hz pulse 
frequencies, the volumes add up to substantial flow ra¬ 
tes. 

The described mechanism can be easily disturbed by 
changes to the dynamics of the colliding columns due to 
dissipation effects, channel morphology and other fac¬ 
tors. In order to optimize pump operation, it is impor¬ 
tant to understand the specifics of expansion, collapse 
and post-collapse phases of the pump cycle as much 
as possible. Advancing a previous study m, we report 
here first successful imaging of single-pulse inertial flow 
captured with a 4 /rs resolution. Consistent with the¬ 
oretical predictions [num, we observe a characteris¬ 
tic N-shape flow in the long arm of the channel. The 
measured flow rate is used to calibrate the full three- 
dimensional (3D) CFD model [T7] and an effective one¬ 
dimensional (ID) model [18] of inertial pumping to ex¬ 
tract the unknown bubble strength. 

Then the frequeney dependenee of inertial pumping 
is studied both experimentally and numerically for the 
first time. A typical pump cycle lasts ^ 50 /is in our 
system. As the pump frequency exceeds 20 kHz, flow 
from a previous pulse is not fully extinguished by the 
time the next vapor bubble is created. Understanding 


temporal interaction of vapor bubbles is of fundamen¬ 
tal importance. We find that at high frequencies the 
overall flow rate continues to increase but pump effi¬ 
ciency (displaced volume per pulse) decreases gradu¬ 
ally. This behavior has to be taken into account when 
designing complex microfiuidic networks driven by in¬ 
ertial pumps. 

2 Methods 

2.1 Experiment 

The microfiuidic system under study consisted of an 
array of 528 U-shape channels 22 /im wide and 17 /im 
tall with a linear density of 600 per inch (see figure [T|). 
Each U-shaped channel had a 33.3 x 15.0 jim^ micro¬ 
heater (a thermal inkjet resistor) located near one end 
of the channel. The channels were fabricated using pho¬ 
toresist (SU8) patterning technique on top of a silicon 
substrate with pre-fabricated microheaters and CMOS 
driving electronics |2T|. Both ends of U-channels were 
fluidically connected to the main fluid reservoir via slots 
etched through the silicon substrate [min]- 

Eabricated devices were filled with a clear test fluid 
of density p = 1.0 g/cm^ and viscosity p = 0.86 ± 0.02 
cP. Viscosity was measured using a Brookfield viscome¬ 
ter at the operating temperature of 45°C. The total 
water content was 70-80%. Microheaters were actuated 
as described earlier in m- The pump frequency could 
be varied from a single pulse to 48 kHz. Elow was ob¬ 
served using a conventional optical microscope with 
coaxial illumination HU using a high intensity xenon 
lamp. Video was recorded with a Vision Research high 
speed camera Phantom vl610 with frame rates between 
250,000 and 1,000,000 frames per second. The exposure 
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Fig. 2 An illustrative high-speed image filmstrip showing tracer movement through a channel. Frame times (shown in top left 
corners) are counted from initiation of the vapor bubble. The vapor bubble is visible in the left leg of the channel in the 4 fis, 
8 fis and 12 fis frames. 


(shutter) time of 450 ns was used. The imaging system 
enabled spatial resolution of 0.8 — 2.5 /im, depending 
on optical magnification used. 

Fluid flows were visualized by 4.0 /im diameter hol¬ 
low polystyrene spheres added to the fluid. Particle im¬ 
ages were processed frame by frame using the camera 
time stamp. An example of particle tracing is shown 
in figure [H Average particle displacement and net flow 
were determined based on multiple pulses for every pump 
frequency tested. Mean displaced volumes and corre¬ 
sponding net flows were calculated by averaging over 
20-30 tracers located at different points in the chan¬ 
nel cross-section. On average, each channel contained 
about one tracer particle at any given time. Thus, sev¬ 
eral dozen pulses had to be processed for every pump 
frequency. Standard deviation of the particles’ displace¬ 
ments across the population was adopted as the dis¬ 
placement error for each frame time. The latter was 
converted to the displaced volume error and flow rate 
error shown by error bars in figures [6] and [71 

High-speed photography and particle image veloci- 
metry (PIV) are powerful tools in microfluidics [231|25]. 
Prior work on bubble-driven pumps was mostly focused 
on the bubble evolution as the vapor-fluid interface pro¬ 
vides good optical contrast Post-collapse 

flow was also visualized using fluorescent imaging m- 
To the best of our knowledge, the present study em¬ 
ploys the highest frame rate to-date to measure flow 
rate over an entire inertial pump cycle. 

2.2 Computational fluid dynamics (CFD) 

CFD modeling was used in combination with theory 
and experiment to gain insight into the behavior of iner¬ 


tial pumps. The first CFD simulations of inertial pump¬ 
ing were reported previously m- In this paper, we use 
the same numerical approach to simulate the physical 
geometry shown in figure [T](a). The CFD code, CFD3, 
is internally developed and is described in detail else¬ 
where [T71[29]. CFD3 is very efficient at modeling mi¬ 
croscale, inkjet flows using the incompressible Navier- 
Stokes equations solved on a staggered, Cartesian grid 
with piecewise linear interface calculation reconstruc¬ 
tion for free surfaces. CFD3 has been used internally 
for modeling of inkjet droplet ejections for many years 
and has been extensively validated with experimental 
data. 

The geometry used in computational analysis is sho¬ 
wn in figure |3l Multiple U-shape loops are in the do- 
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Fig. 3 Three-dimensional isometric drawing of the device 
used by the CFD model. Boundary conditions for CFD sim¬ 
ulations are set at regions indicated by letters A, B, and C. 
The black rectangle is the pump. 
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Fig. 4 A typical CFD run at pump frequency of 36 kHz. 
Time between pulses is 27.8 fis. Main panel: the overall time 
series. Inset: displaced volume as a function of pulse number. 
The volume per pulse converges to a quasi-stationary limit of 
7.2 pL/pulse. 

main and the nozzles (indicated by A in figure [3|) are 
open to the atmosphere which is modeled using a gas 
at atmospheric pressure. The center loop contains a mi¬ 
croheater which is used to simulate boiling through a 
process described in Ref. [171129]. The bottom of the 
domain (indicated by C in figure [3j) is modeled using a 
pressure boundary condition. This pressure boundary 
condition is atmospheric pressure minus a small back 
pressure, which is approximately 1000 Pa in this case. 
All other boundaries have solid walls (indicated by B 
in figure [3]). 

The experimental pumping event is simulated by 
a microheater boiling event and subsequent tracking 
of the integrated amount of flow through the central 
U loop. The fluid in the simulation has viscosity /i = 
0.86 cP, a density p = 1.0 g/cm^, a contact angle 45°, 
and surface tension a = 50 dyn/cm. The amount of 
pressure-volume work performed by the boiling event 
is determined by an internal variable known as critical 
volume or Rcr- Critical volume is an adjustable param¬ 
eter that defines the pressure in the vapor bubble, and 
subsequently, the amount of work done by the vapor 
bubble. It is defined in relative units referencing a re¬ 
sistor operating with a standard input pulse on a ref¬ 
erence fluid. The value of vary from 1 if the 

fluid boils differently than the reference ffuid or a non¬ 
standard input pulse is used. 

In contrast to the previous work [T7] , here we simu¬ 
late multiple successive pump cycles over long continu¬ 
ous runs. Each subsequent bubble is injected into rem¬ 
nant ffow from the previous cycle. After several pulses 
the system reaches a quasi-stationary state with dis¬ 
placed volume being independent of pulse number. This 
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Fig. 5 Simplified ID schematic of an inertial pump. The 
small black rectangle marks the location of a microheater xq. 
xq is also the initial positions for both vapor-fluid interfaces 
at the start of bubble expansion, po is the external pressure 
equal in both reservoirs. The residual vapor pressure pvr = 
0.3 atm is assumed to be independent of xq. pi and p 2 are 
the dynamic pressures at the channel-reservoir interfaces and 
are inflow-outflow asymmetric. 

is a more accurate way of assessing interaction between 
pulses. In m, pump frequency effects were estimated 
by truncating single pulse ffows at appropriate times, 
the single pulse always starting from a state of rest. 

A typical CFD3 pumping result for multiple pulses 
is shown in figure 01 The cumulative amount of ffuid 
volume passing through the central U-loop is measured 
and shown as a function of time in the main panel. The 
next step is to calculate displaced volume per pulse as 
shown in the inset graph. Typically, within a few pulses, 
the displaced volume approaches a quasi-stationary va¬ 
lue, which is 7.2 pL/pulse in this example. 

2.3 The one-dimensional model 

A simplified one-dimensional model of inertial pumping 
was developed previously m- The basic geometry is 
shown in figure [H A pump state is completely described 
by positions of two vapor-ffuid interfaces xi (t) and X 2 (t) 
and their temporal derivatives xi{t) and X 2 {t)- In terms 
of dimensionless time r and interface coordinates ^ 1 ^ 2 , 
cf. table [U the dynamic equations for the left and right 
ffuidic columns read 

+ = 1, (1) 

(1 - 6) ^2 - 1^2 H{-Q + /3 (1 - 6) = - 1 . (2) 

Here the prime denotes a derivative with respect to r, /3 
is effective friction, and H{z) is a Heaviside step func¬ 
tion: H{z > 0) = 1, and H{z < 0) = 0. The pres¬ 
sures at the channel-reservoir interfaces are given by 
Pi ,2 = Po — ^pxi 2 H{±xi^ 2 )- The interface pressures 
are equal to the external pressure po during expansion, 
when the ffuid is ejected into reservoirs. During col¬ 
lapse, the pressure is reduced by a Bernoulli term m- 
Using CFD3, we have verified that the ffow-asymmetric 
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Variable / parameter 

Formula 

Values 

Time 

- t /Po-Pvr 

^ V p 

0-5 

Interface position 

II 

67 

(0.0 - 1.0) 

Interface velocity 

j.. - ^ - P' 

Ot 

0-1 

Microheater location 

0 

II 

(0.0 - 1.0) 

Bubble strength 

a- ^ 

pAL V Po-Pvr 

0.1 - 1.0 

Effective friction 

0 _ npL / p 

^ pA V Po-Pvr 

1-10 


Table 1 Dimensionless variables and parameters used in the 
effective ID model of inertial pumping US]. Here t is physical 
time and xi ,2 are the physical locations of vapor-liquid inter¬ 
faces. xq is the physical location of the microheater center, L 
is the channel length, A is the channel cross-sectional area, p 
and pi are the fluid density and viscosity, po = 10^ Pa is the 
atmospheric pressure, and pvr = 0.3 • 10® Pa is the residual 
pressure of the vapor bubble. Quantity qo is the initial liq¬ 
uid momentum that is approximately given by the product 
of vapor pressure, channel cross-section and duration of the 
high-pressure phase, k is the geometrical loss factor that de¬ 
pends on the shape of channel cross-section. = 29.34 for a 
22/17 rectangle (the present case) |30| . The last column lists 
values typical for the experimental systems studied in this 
work. 

dynamic pressure is more representative of the inertial 
pump physics than the constant boundary conditions 
also considered in Ref. m- After the two columns col¬ 
lide at point at time Tc, the entire channel fluid moves 
as a whole with a constant mass. Fluid motion in this 
post-collapse phase is described by an equation 

r + ie'"sgn(n+/3r = o, (s) 

with ^(tc) = and ^'(tc) = r/c- The differential equa¬ 
tions (HD-® can be solved by direct integration in time 
domain using standard numerical methods. 

The initial conditions for the expansion-collapse ph¬ 
ase are 

6(0) = 6(0) = ^0, (4) 

6(0) = -| + , 6(0) = + m ■ (5) 

where a is the dimensionless bubble strength, is the 
microheater location, and r]o is the velocity prior to 
bubble generation. Dimensionless model parameters are 
defined in table [H Note that the starting point of bub¬ 
ble expansion (which coincides with the microheater 
location) can be any value on the interval (0,1), only 
excluding end points 0 and 1. Mathematically, the dy¬ 
namic model is well defined even in the limit ^ 0- 
A small fluid mass results in a large deceleration that 
immediately turns around the flow, so that the singular 
point ^ = 0 is never reached. Clearly, in real devices 
is limited by finite sizes of the microheater (20 — 50 jam 


in length) and of the microchannel (L = 200 — 2000 
jam in length). Thus, in practical calculations the start¬ 
ing point ^0 = xq/L typically lies within the interval 
0.01 < Co < 0.99. 

A non-zero starting velocity r]o is a novel feature of 
the ID theory compared with the single-pulse case stud¬ 
ied in Ref. [18]. If a new pulse is initiated before flow 
from the previous pulse is completely extinguished, the 
new vapor bubble will be initiated in a moving fluid. 
This changes the expansion-collapse dynamics of both 
fluidic columns and may affect overall net flow. The 
pre-bubble velocity r]o has a one-to-one correspondence 
with the pulse frequency /. In practical calculations, 
770 is set to some value, then expansion-collapse flow is 
solved to obtain post-collapse velocity rjci and finally 
post-collapse flow is solved to find the moment when 
flow slows down from pc to 770 - The procedure deter¬ 
mines a total time after which the system returns to 
its original state. An inverse of this time, expressed in 
physical units, is the pump frequency that can be di¬ 
rectly compared with experiment. 

We close this section by comparing the benefits of 
numerical and theoretical methods. CFD provides the 
most accurate treatment of an inertial flow problem, 
and is typically the first method to try on a new sys¬ 
tem. It takes into account flow features in reservoirs, 
finite microheater dimensions, realistic vapor pressures 
and other elements of a real architecture. However, sim¬ 
ulations take a long time: runs reported in this work 
take up to several days to complete. In contrast, the 
ID model equations can be solved in a few seconds, but 
they do not capture the system elements listed above. 
The ID model always requires calibration with both 
experiment and CFD to make sure the model parame¬ 
ters are realistic and in correspondence with practical 
devices. Once the key parameters are set (for exam¬ 
ple the unknown bubble strength a), the ID model be¬ 
comes an excellent tool for exploring the vast space of 
microfluidic architectures. We believe that at this stage 
of inertial pump development, it is necessary to carry 
both approaches as they are complementary. 

3 Results 

First, we characterized individual pumping events at 
low pump frequencies below 10 Hz. At such frequen¬ 
cies, successive pulses do not overlap in time. Volume 
displaced through the long channel arm versus time is 
plotted in figure [6l Experimental data points are shown 
by circles. Note that the curve has a characteristic N- 
shape. Once the vapor bubble is created, the fluid is 
“kicked” through the channel by high vapor pressure. 
Once the vapor cools off and the pressure drops (within 
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Fig. 6 Single-pulse flow curve. Circles with error bars are 
experimental data. The solid line is CFD flow for Vcr = 1.1. 
The dashed line is the best fit by the ID model with param¬ 
eters fo = 0.158, a = 0.276. The experimental error bars 
represent standard deviations of the displacement of about 
30 tracer particles located at random points of the channel 
cross-section. 


a couple of /is), the flow starts to decelerate under the 
combined action of friction forces and a negative pres¬ 
sure difference. At around 8 /is, the long arm flow re¬ 
verses direction and the bubble collapses at about 12 
/is. During collapse the long arm has a smaller mechan¬ 
ical momentum than the short arm (for the reasons ex¬ 
plained in Refs. nain) and the flow reverses direction 
again, now from negative to positive. This secondary, 
post-collapse flow lasts for another 40-50 /is until it is 
fully extinguished by friction forces and the entire fluid 
comes to rest. Note that the N-shape flow was predicted 
theoretically iniiiH], but it has been observed experi¬ 
mentally for the first time here. 

Results from the CFD simulations are presented in 
figure [6l Comparisons to experimental data were used 
to determine the critical volume Vcr = 1-1. The fluid 
parameters used in the simulations were described in 
section 12.21 The CFD simulations match the experi¬ 
ment quite well during the initial phase until the bubble 
collapses at about 12 /is. Then there is a small devia¬ 
tion between the experiment and numerical model. The 
small hump in displaced volume in the CFD simulation 
from 14 to 18 /is is due to bubble rebound. This re¬ 
bound is seen in our experiment (and it was observed 
in the past miisi]) but its effect on displaced volume 
is masked by the volume measurement error. 

The best ID model fit of single-pulse flow is shown 
in figure [6] by a dashed line. The overall agreement is 
good, except that the total displaced volume comes out 
slightly larger than the measured value (8.2 pL vs 7.8 



Fig. 7 Experimental flow rates at different frequencies su¬ 
perimposed on predictions of CFD (for Vcr = 1-1) and the 
ID model with parameters extracted from single-pulse fits. 


pL) but still within experimental error bars. In fitting, 
the channel length (403 /im), width (22 /im), height 
(17 jim) as well as fluid density p = 1.0 g/cm^ and 
viscosity p = 0.86 cP were set to experimental values. 
The microheater location = xojL and effective bub¬ 
ble strength a were left as adjustable parameters. The 
optimal values have been found to be = 0.158 and 
a = 0.276. The a value falls right within the estimate 
given in table [H = 0.158 is about 10% larger than 
the “theoretical value” of 0.140. The latter is obtained 
if the xq is associated with the geometrical center of 
the pump microheater, which is 56.5 pm away from the 
reservoir, cf. figure [H The remaining differences can be 
attributed to a finite size of the microheater, presence 
of open nozzles and particle filters that are not captured 
within the ID model. 

The flow rate dependence on pump frequency is now 
discussed. Quasi-steady flow rates were measured at 
/ = 12, 18, 24, 32, 36, 42, and 48 kHz. Experimental 
data points are shown in figure 0 by circles. Since the 
duration of one pulse is about 50 /is, successive pulses 
do not overlap in time until / exceeds 20 kHz. Thus, 
the flow rate is expected to be linear in / for / < 20 
kHz. At higher frequency, the flows begin to interfere 
with each other. It is of fundamental interest to know 
the sign of the effect. Experimental data suggest that 
the interaction is negative, i.e. the flow rate grows with 
the pump frequency sub-linearly. 

The sublinear behavior at high frequencies is repro¬ 
duced by both CED and ID models as shown in flgure0 
In both cases, model parameters were taken from the 
single pulse fits described above. In the CED case, the 
run time was long enough to reach the quasi-stationary 
flow, as explained earlier. There is excellent agreement 
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Fig. 8 Steady-state volume per pulse compared with flrst- 
pulse volume for different pump frequencies. The lines are 
parabolic fits to guide the eye. 

between CFD and experiment, within experimental er¬ 
ror bars. The ID model slightly overestimates the flow 
at high frequencies. We attribute this to be a conse¬ 
quence of the single-pulse fit that produced a displaced 
volume slightly higher than experiment, cf. figure [6l 
That may have led to a slightly larger effective bubble 
strength and less interaction between successive pulses 
as a result. Overall agreement between our models and 
experiment is good. 

Physically, the reduction of flow rate is caused by 
the new vapor bubble disrupting post-collapse flow gen¬ 
erated by the previous pulse. For a well-designed pump, 
post-collapse flow provides the bigger portion of the 
overall displaced volume, and the new bubble cuts short 
the contribution of this secondary phase to the previ¬ 
ous pulse. As a result, the volume per pulse gradually 
goes down. It is instructive to compare the effect of 
such short-cutting on both first-pulse and steady state 
flows, which is done in figure [H The first pulse occurs in 
a stagnant flow and the corresponding displacement is 
defined by the moment when the single-pulse flow curve 
of figure [6] is terminated by the second pulse. However, 
subsequent pulses take place in nonzero flow, which in¬ 
creases the overall displacement for the same time in¬ 
terval between pulses. As a result, some volume lost at 
first pulse is recovered as the flow approached steady 
state. That is why the steady state curve shown in fig¬ 
ure [8] is more flat that the first-pulse curve. 

4 Summary 

Due to its small size and complete integrability, the in¬ 
ertial pump may become an essential element of future 
fluidic integrated circuits. It is important to understand 


various aspects of this new technology. In this paper, 
single-pulse dynamics has been imaged by a high-speed 
camera for the first time. This enabled calibration of the 
full three-dimensional CFD model and an effective one¬ 
dimensional model with respect to the unknown bubble 
strength. Then, variation of time-average flow rate with 
pump frequency was investigated, also for the first time. 
Both experiment and models have shown sublinear be¬ 
havior at high frequencies. Such behavior is caused by 
temporal interaction of successive pump pulses, when 
each new pulse disrupts the post-collapse flow gener¬ 
ated by the previous pulse. As a result, the overall flow 
rate still grows at high frequencies but the displaced vol¬ 
ume per pulse goes down. From this qualitative picture, 
the sublinear dependence of the total flow rate should 
be more prominent in systems with longer post-collapse 
flows, e.g. in microchannels with larger cross-sectional 
areas than those studied in this work. 
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